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A detailed zero and transverse field (ZF&TF) muon spin rotation (/xSR) investigation of mag- 
netism and the magneto-electric coupling in Cu2 0Se03 is reported. An internal magnetic field 
Bi-nt{T = 0) = 85.37(25) mT was found, in agreement with a ferrimagnetic state below Tc — 57.0(1) 
K. The temperature dependence of the magnetic order parameter is well described by the relation 
-Bint = -B(0)(1 — {T /Tc)"^)^ with an effective exponent /3 ~ 0.39(1) which is close to the critical ex- 
ponent /? ~ 1/3 for a three dimensional (3D) magnetic system. Just above Tc the muon relaxation 
rate follows the power low A(r) oc {T /Tc — 1)^" with Cb — 1.06(9), which is characteristic for 3D fer- 
romagnets. Measurements of i3int(r) with and without an applied electrostatic field E = 1.66 x 10'' 
V/m suggest a possible electric field effect of magnitude ABv = -Bv(0V)-i3v(500V) = -0.4(4) mT. 
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I. INTRODUCTION 



Much attention has been directed toward muhiferroic 
and magnetoelectric materials in recent yearsii"— The 
coupling between magnetic and electric parameters in- 
creases the degrees of freedom of the ordered ground 
state, making these materials good candidates for the 
study of new phenomena in highly correlated electronic 
systems. Strong magnetoelectric coupling is rather rare 
in the solid state, since usual microscopic mechanisms for 
magnetic and electric polarization are mutually exclu- 
sive. Magnetism requires strong exchange interactions 
related to a strong hybridization of the transition ion 
electrons leading to conductivity. Conductivity, on the 
other hand, is inconsistent with the presence of an elec- 
tric polarization in a samplei^ It is therefore of particu- 
lar importance to unravel the mechanisms behind mag- 
netoelectric coupling. A number of atomic mechanisms 
have been proposed in order to explain the magnetoelec- 
tric couplingjii^^— Considering spatial and time inversion 
symmetry for the magnetic (M) and the electric (P) po- 
larization, it was concluded that linear magnetoelectric 
coupling is only possible when both vectors vary in space 
and time4 On the other hand, the importance of frus- 
tration effects in magnetoelectrics for the static polar- 
izations of P and M was stressed for nonlinear coupling 
mechanisms 1^ The presence of large magnetic and elec- 
tric polarizations is an important condition for strong 
magnetoelectric coupling, making ferro- or ferrimagnetic 
materials favorable candidates li 

The ferrimagnetic magnetoelectric compound 
Cu2 0Se03 was recently discovered and single 
crystals were successfully grown soon afteriii The 
compound is piezoelectric and undergoes a ferrimag- 
netic transition below 60 K, exhibiting magnetoelectric 
coupling as revealed by magneto-capacitance studies 
on a polycrystalline sample.— An abrupt change of the 
dielectric constant below the ferrimagnetic transition 
was later confirmed by infrared studies J^i^ At present 



the nature of the magnetoelectric coupling is unknown, 
since the most common mechanisms, involving magne- 
tostriction and piezoelectric effects via lattice distortions 
are excluded. Neither X-ray diffraction (XRD)^ nor 
nuclear magnetic resonance (NMR)ii studies revealed 
any lattice anomaly below the Neel temperature. 

The positive muon is a very sensitive microscopic probe 
for studying magnetic properties in zero as well as in an 
applied electric field Following the pioneering works of 
Eschchenko et alJ^ and Lewtas et alt^ we implemented 
a setup with alternating electric fields, and performed 
a muon spin rotation/relaxation (/xSR) investigation of 
magnetism and magnetoelectric coupling in Cu2 0Se03. 
The temperature dependence of the internal magnetic 
field -Bint was investigated below Tc = 57.0(1) K, and 
the relaxation rate was studied above T^. Zero field /xSR 
measurements of the internal field distribution with and 
without an applied electric field E = 1.66 x 10^ V/m indi- 
cate a small electric field effect on the internal magnetic 
field: ABy = By(OV) - By (500 V) = -0.4(4) mT. 

The paper is organized as follows: The sample prepa- 
ration and the details of experimental setup are described 
in Sec. [Til In Sec. |lll]we describe the model used for 
the analysis of the /iSR data and the relation of the mea- 
sured /iSR spectrum to the lattice and magnetic struc- 
ture of the sample. In Sec. IIVI we present and discuss 
the obtained results. Conclusion is given in Sec. |Vl 



II. EXPERIMENTAL DETAILS 

A high quality single crystal of Cu20Se03 of ap- 
proximate size 7x8x3 mm^ was prepared in a manner 
described elsewherej^i The zero field (ZF) and trans- 
verse field (TF) /iSR experiments were performed at the 
7rE3 beam line at the Paul Scherrer Institute (Villigen, 
Switzerland). The crystal structure is cubic with sym- 
metry (P2i3)i^iiS The sample was aligned with its (100) 
direction parallel to the incident muon beam. The spin 
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vector of the muon was oriented approximately with an 
angle of 60 degrees with respect to the momentum vec- 
tor. The asymmetry time spectra were monitored in the 
"Forward", "Backward", "Up", and "Down" (FBUD) 
positron detectors Typical statistics was 40 to 50x10® 
positron events in the FBUD detectors for a spectrum 
with a 5 /is time window. The switched electric field E 
was applied along the (100) direction of the 3 mm thick 
crystal. The crystal was mounted between two Cu elec- 
trodes: A thin (50/^m) Cu metal foil was used as the posi- 
tive electrode, and the negative electrode was soldered to 
the sample holder and electrically connected to ground. 
The applied voltage was switched at a rate of 100 Hz for 
two reasons: (i) to avoid accumulation of muon created 
charge in the vicinity of the sample-to-electrode contact 
surface, which might offset and even cancel the applied 
field E, and (ii) to provide a well-defined consecutive ref- 
erence in order to reduce artefacts related to any slight 
variation in temperature or applied magnetic field with 
time. All the positron events registered by the FBUD 
detectors were stored alternatively in the first block of 
four histograms when the electric field was off, and the 
second block of four histograms when the electric field 
was on. Measurements were performed for two different 
electric fields: E = 500/3 V/mm and E = 800/3 V/mm 
(i.e. 500 or 800 V applied on 3 mm thick sample). 
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FIG. 1: (Color online) /iSR time spectra of Cu20Se03 at 10 
and 58 K in ZF with (0) and without (o) applied electrical 
field E = 500/3 V/mm. Solid lines are fits to the data using 
Eqs. m and ®. 



III. ANALYSIS AND MODELS 

ZF /iSR allows to determine internal magnetic fields 
at the position in the lattice where the muons stop. For 
a polycrystalline sample with a given static magnetic 
field B at the stopping site of the muon, the muon de- 
polarization functioni^ may be expressed as: Gsit) = 
\ + l cosin ^Bt), where 7^ = 27r x 135.53 MHz/T is the 
gyro- magnetic ratio of muon. Gsit) consists of two parts: 
a constant fraction of 1/3, and a fraction 2/3 oscillating 
with the frequency uj = 7pi?. For a given magnetic field 
distribution P{B), the muon depolarization function at 
the muon site isi^i 
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Gp(t) - - + - P{B) cosi-f^Bt) dB. (1) 

The function P{B) contains information on the magnetic 
structure and the spatial magnetic field distribution of 
the sample, but also the effect of a static electric field E 
on the local magnetic fields in the sample. 

Due to an additional isotropic dynamic muon relax- 
ation and the presence of a Cu background signal, the ZF 
depolarization function for Cu2 0Se03 is best described 
by the following equation: 



A(t)-As-Q + ^^ P{B) cos{j^Bt + ^) dB 



+ABg ■ GKrio'Bgt) 



(2) 



Here is proportional to the fraction of muons stopping 
in the sample, and Asg is proportional to the fraction of 




0.00 -a 



B(G) 



FIG. 2: (Color online) Fourier transform of the oscillating 
part of the asymmetry spectra shown in Fig. [1] (black circles: 
E = 0; red diamonds: E — 500/3 V/mm). The solid lines are 
the FT of the corresponding solid curves in Fig. [T] 



muons stopping in the Cu sample holder. The parameters 
Xsg and asg describe the temperature independent muon 
depolarization in Cu, while A corresponds to the dynamic 
muon relaxation in the sample. GxTis^Bgi) denotes the 
Gaussian Kubo-Toyabe depolarization function.— Note 
that the phase = for zero field. As for a powder sam- 
ple, we observe also for the present single crystal sample 
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a static fraction of 1/3 for the depolarization function 
and an oscillating fraction of 2/3J^ In the case of a sin- 
gle crystal sample, though, these fractions are due to the 
formation of magnetic domains with a random spatial 
distribution which in fact corresponds to the situation 
for a polycrystalline sample. The best fit of P{B) to the 
experimental data with a minimal set of parameters was 
obtained with a sum of two Gaussians: 



■ exp 



B - B, 



(3) 



Here Fg and Fj arc the fractions of the two Gaussians 
with mean fields Bi and B2 and standard deviations 
f 1/7^1 Analysis of the data measured with 

highest statistics (100 x 10^ positron events) at 10 K 
with Eqs. ^ and (O [see Figs. (P and ©] yields 
in Fl = 0.18(2), F| = 0.82(2), and B2/B1 = 1.07(2) 



[note that F^ + Fg = 1]. The fit were performed by 
keeping all the sample parameters (^5, cri, (T2, A, i^J, 
Fg, and B2/B1) and the background parameters (^sg, 
(Tsg, and Xbq) the same, while i3i(0V) and i?i(500V) 
(data without and with the applied electric field) were 
free parameters. The temperature independent sample 
asymmetry As = 0.144 and the background parameters 
ABg = 0.05, (JBg = 0.14(2) /is-\ and Xsg = 0.11(2) 
/is~^ were determined from the global fit of the whole 
temperature dependence, while the total initial asymme- 
try As-\-Abg = 0.194 was determined above T^. Further 
more, all the background parameters and the following 
sample parameters: Fg, Fj, B2/B1, were kept as tem- 
perature independent. 

Figure [T] shows the asymmetry time spectra at 10 K 
and 58 K (below and above the magnetic transition) with 
and without electrostatic field E. The corresponding fits 
using Eqs. 1^ and are represented by the lines. Fig- 
ure [2] shows the Fourier transform (FT) amplitudes of 
the oscillating part of the /xSR spectra and the fitted 
curves shown in Fig. [T] For a small relaxation rate A 
these FT amplitudes represent the magnetic field distri- 
bution P{B) given by Eq. Note that P{B) consisting 
of two Gaussians [Eq. ([3])] describes the basic features of 
the /xSR spectrum quite well, and that the FT amplitudes 
for F = and E = 500/3 V/mm are almost identical. 

Before presenting the experimental results, we describe 
below the relation between the measured /xSR spectra 
and the magnetic structure of Cu20Se03. The crystal 
symmetry of Cu2 0Se03 is cubic with a lattice constant 
a = 8.9IIIA and spatial group symmetry P2i3i^ The lat- 
tice structure is thus quite complex with 32 oxygen ions in 
the unit cell, suggesting quite a large number of possible 
muon stopping sites. Generally, the positive muon stops 
at a high symmetry interstitial site of the lattice, close to 
negatively charged ions (in this case O^^). In some cases 
the muon may form weak bonds with oxygeufi^ The func- 
tion P{B) describes the muon site-weighted distribution 
of internal fields in the sample. 

In order to find the muon stopping sites, an analysis of 
the electrostatic potential was performed. The potential 
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FIG. 3: (Color online) Magnetic structure and dipolar fields 
in Cu20Se03 (a). The black spherical surfaces represent a 
field strength of 20 T (i.e. Cu'^+ ions). The arrows indicate 
the magnetic moments of Cu'^"''. Blue surfaces represent in- 
ternal fields of 10 mT (for 0.5/xs per Cu ion). The red spots 
A, B, C, D, and E are the muon stopping sites. Panel (b) 
shows the same as (a) but in the (111) direction. 



at position r within the lattice unit cell was approximated 
using a spherical point charge-like model (in dimension- 
less units): 



V(r)=E 
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FIG. 4: (Color online) Probability distribution of dipolar 
fields Pdip(B) with equal weights for sites A, B, C, D, and 
E (•). The contributions of the individual sites A, B, C, D, 
and E to Pdip(-B) are represented by the symbols, V, A, o, (), 
and □, respectively. The width of each curve is proportional 
to the magnetic field gradient at the corresponding site. The 
dotted line shows the Pdip(-B) with an additional contact field 
of 70 mT at site E. 



The first term is the Coulomb potential, while the sec- 
ond term is the exchange correlation potential, which is 
usually assumed to be proportional to the local charge 
density (the expression exp(— represents an av- 



erage radial wave function" 
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Here, Rj = |r — rJ, and 



the sum is taken over ion coordinates rj within a clus- 
ter of 4 X 4 X 4 unit cells. The charges qj and the ionic 
radii Rj are: +2 and 0.71 A, +4 and 0.42 A, -2 and 
1.3 A for Cu, Se, and O ions, respectively (here the ele- 
mentary charge is unity). The adjustable parameter V^^ 
was chosen as ~ -f-lO. We found that the coordinates of 
the potential minima do not appreciably depend on V^^ 
over a broad range of values. The potential V{r) has five 
magnetically non-equivalent minima with nearly equal 
depth at the following sites: A = (0.215,0.700,0.970), B 
= (0.035, 0.720, 0.805), C = (0.195, 0.555, 0.795), D = 
(0.275, 0.295, 0.460), and E = (0.635, 0.550, 0.525). The 
sites A, B, C, D, and E arc indicated as red spots in Fig. 
[31 In addition, there arc four local minima with higher 
energy and lower probability to be occupied. 

The muon probes the vector sum of the internal (dipo- 
lar) magnetic field and the contact field at a particular 
lattice cite. The dipolar magnetic field B(r) at position 
r within the lattice unit cell is calculated as follows)^ 



i,/3 




/3 



gap 



(5) 



Here Ri = r — r^, a and (3 denote the vector components 



X, y, and z, is the position of i-th magnetic ion in the 
unit cell, and is the corresponding dipolar moment. 
The summation is taken over a sufficiently large Lorcntz 
sphere of radius R^. Beyond the Lorcntz sphere, the in- 
tegration is carried out over the domain volume. The 
contribution to the internal magnetic field from this in- 
tegral is B' = 47r//o(3 - Nd)]VI^ where //qM ~ 66 mT 
is the domain magnetization, and the demagnetiza- 
tion tensor determined by the geometry of the domain 
and the magnetic anisotropy. For the calculation of the 
magnetization the lattice constant a = 8.91113 A and 
the magnetic moment of 0.5/iB per Cu^+ ion were taken 
from Ref. For a magnetically isotropic spherical do- 
main Nd = i.— The field B' and the stray fields due 
to the neighboring domains average statistically to zero. 
They only give rise to an additional broadening of P{B), 
which is smaller than the width of the narrow compo- 
nent of P{B) (see Fig. The magnetic structure of 
CU2OSCO3 and the spatial magnetic field distribution 
calculated with Eq. ([5]) is shown in Fig. [31 The prob- 
ability field distributions for the magnetic structure with 
equal weights for the muon sites A, B, C, D, and E are 
shown in Fig. 31 The calculations were performed with 
Gaussian sampling around the points A, B, C, D, and E, 
with a standard deviation ctl — 0.23 A. Thus, the widths 
of the curves in Fig. [H are proportional to the magnetic 
field gradients at the corresponding sites. The total field 
distribution from all sites A, B, C, and D has a broad 
Gaussian-like shape centered at approximately 100 mT. 
This broad distribution agrees quite well with the experi- 
mental distribution (see Fig. [3]). The narrow peak calcu- 
lated for site E is located at about 15 mT, in contrast to 
the experimentally observed peak at around 85 mT. This 
discrepancy may be explained by assuming an additional 
contact field of approximately 70 mT at the muon site E, 
resulting in a peak position of 85 mT (the total field is 
the vector sum of the dipolar field Bdip and the contact 
field Be). Note that the ratio of the broad and the nar- 
row signal intensities is about 4, in good agreement with 
the ratio Fg/Fg ^ 4.5(5) obtained from Eq. above. 



IV. RESULTS AND DISCUSSION 

The temperature dependence of the parameters cti , (T2 , 
and A as obtained from the data analysis by means of Eqs. 
([2]) and ^ are displayed in Fig. [5l The parameters (Ti 
and (72 decrease with increasing temperature and drop 
to zero at the Curie temperature Tc = 57.0(1) K of the 
ferrimagnctic transition of Cu20Se03. The longitudinal 
relaxation rate A is a measure of the internal magnetic 
field dynamics and can be expressed by the field-field cor- 
relation function A = 7^ < B±{t)B±{0) > dtM. Here 
B±^ is the field component perpendicular to the muon- 
spin direction at the muon site. Brackets denote statis- 
tical averages. In the paramagnetic state close to the 
relaxation rate A is a measure of the spatial correlation 
length ^ of the magnetic order. In the critical state the 
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FIG. 5: (Color online) Transverse relaxation rates ai and 
(72 as well as longitudinal relaxation rate A of single-crystal 
Cu20Se03 as a function of temperature. The black solid 
line represents a fit of the data to the power law A(T) = 
A{T/Tc - 1)"" (see text). 



following relations hold: A oc ^'^/^ for fcrromagncts and 
A cx ^^/^ for antiferromagnetsi^"— Note that these re- 
lations are only strictly valid in the critical regime very 
close to Tc- However, it was found empirically, that these 
relations describe experimental data rather well for EuO, 
EuS, RbMnFg up to ~ LST^,^^ Above Tc the correla- 
tion length ^ follows the power low ^ oc {T/Tc — 1)^'', 
where the effective critical exponent D ~ 0.70 for a 3D 
Heiscnberg magnet j^"— A fit of A above Tc to the equa- 
tion A(T) = A(T/Tc - l)-"^ results in A = 0.0030(5) 
fis~^ and the effective dynamic exponent2S J) — 1.06(9) 
(see the black line in Fig. [5]). Thus, for i> ~ 0.70, we find 
with u/ij = 1.5 that A oc in fair agreement with the 
critical behaviour of a 3D ferromagnet (A oc S,^^'^)- 

Next we discuss the fiSK experiments with an electric 
field applied. The temperature dependence of the mean 
internal magnetic field Bi with or without applied elec- 
tric field E is shown in Fig. [DJa). The temperature de- 
pendence of the internal magnetic field (i.e. the magnetic 
order parameter) decreases with increasing temperature 
and vanishes at Tc. Analyzing the data with the power 
lawi^ 



B,=B,{Q)-il-{T/TX)^ 



(6) 



yields Tc = 57.0(1) K, a = 2.00(9), P = 0.39(1), and 
Bi{0) ^ 85.37(25) and 85.57(25) mT for S = and 
500/3 V/mm, respectively. The value of the effective 
critical expoment^ /3 lies quite close to the critical expo- 
nent /3 ~ 1/3 expected for a 3D magnetic system.— For 
T the temperature dependence of the internal field 
Bi{0) — Bi{T) oc T" is determined by the excitation of 
the ground state magnetic order i22i2^ Figure HJb) shows 
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FIG. 6: (Color online) (a) Temperature dependence of the 
mean internal Bi field in single crystal Cu20Se03 for zero 
and applied electrostatic field E. The solid lines are power 
law fits to the data with Eq. The insert shows Bi vs 

(1 - {T/Tcf) on a log-log scale.(b) Electric field shift ABy = 
B\ (OV) — B\ (500V) as a function of temperature. The solid 
red line corresponds to the mean value of ABv — —0.4(4) 
mT below 50 K. (c) Electric field shift ABy = BziOV) - 
i?2(800V) as a function of magnetic field measured in the TF 
fiSK experiment at lOK. 



the difference ABy = Bi(0V)-Bi(500V) as a function of 
temperature. Obviously, the data points scatter substan- 
tially. However, a close examination shows that the data 
points below 50 K are systematically shifted to negative 
values, except of a point at 10 K [the blue open circle 
in Fig. [6lb)] which is the first measured point and it 
was recorded in the test phase of the experiment. Thus, 
we cannot exclude additional systematic errors related 
to instrument set up for the first point. The statistical 
average of all values of ABy is zero within error. By 
excluding the data point at 10 K the average value ABy 
below 50 K was found to be ~ —0.4(4) mT. These exper- 
iments suggest the presence of a small electric field effect 
on the magnetic state of Cu20Se03. 

Additional fiSK experiments were performed with an 
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applied transverse magnetic field (TF). The ^SR spec- 
tra were found to be well described by the asymmetry 
function: 

ATF{t) :=isexp (^"^^It^^ cos(7^B2i + (^)+ 

+ABgexp (^-\^lgt^^ cos(7^SBgi + 0). (7) 

The parameter As{ABg) corresponds to As{ABg) in Eq. 

while <T2{(JBg) and B2{BBg) are the relaxation rate 
and the mean field of the sample (Cu background), re- 
spectively. The values of 172 and (T2 were found to be 
approximately the same, and the Cu relaxation rate 
(TBg — 0.25 /is~^ is small. ^2 is slightly larger than the 
applied field, while is close to the applied field. In 
the TF /iSR experiments the electric field amplitude was 
increased to 800/3 V/mm. Note that the signal from the 
sample is well described by a single Gaussian, and that 
the narrow signal in Eq. ([3]) broadens. The field depen- 
dence of AB'y = i32(0V) - i?2(800V) at T = 10 K is 
shown in Fig. [DJc). No effect of the electric field E on 
the internal magnetic field was found within the precision 
of the experiment, although a strong field dependence of 
the magneto-capacitance was reported.— The TF ^SR ex- 
periment is less precise than the ZF experiment, since the 
narrow signal that mainly determines the errors becomes 
broader by applying a magnetic field. 

V. CONCLUSION 

In conclusion, we performed a detailed investiga- 
tion of magnetism and the magneto-electric effect in 



Cu20Se03 by ZF and TF ^SR. An internal magnetic 
field Bi„t(T = 0) = 85.37(25) mT was detected below 
Tc = 57.0(1) K, consistent with a ferrimagnetic state.— 
The effective critical exponent for the temperature 
dependence of i?int was found to be ^ ~ 0.39(1), in fair 
agreement with the critical exponent /3 ~ 1/3 expected 
for 3D magnetic systems. The magnetic order parameter 
B{0) — B{T) (X was found to exhibit a quadratic 
temperature dependence for T — > 0. The temperature 
dependence of the muon relaxation rate above Tc is 
well described by the relation A cx {T/T^ — 1)^" with 
Cb ~ 1.06(9), suggesting A (x ^"^/^ (^ is the magnetic 
correlation length) which is characteristic for 3D 
ferromagnets .^^1^^ The ZF /iSR measurements of the 
microscopic internal field distribution with and without 
applied electric field E — 500/3 V/mm indicate a 
small electric field effect on the internal magnetic field: 
ABv = Bi(OV) - Bi(500V) = -0.4(4) mT. The strong 
muon relaxation sets a limit on the precision of detecting 
a magneto-electric effect. To improve the precision of 
the ^SR experiment substantially higher statistics would 
be needed. 
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